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The formation of a long-range Jahn- Teller distortion in undoped manganites is studied within 
the dynamical mean-field theory taking both electron-phonon and electron-electron interactions 
into account. We find that the observed insulating behavior at temperatures above the Jahn- Teller 
transition and the bond-length changes below it can most naturally by explained by assuming a 
strong electronic repulsion and a rather weak electron-phonon coupling. 
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The physics of manganese oxides Rei-^A^MnOa (Re 
is a rare earth metal ion such as La or Pr, A a divalent 
alkali such as Sr or Ca) is characterized by strong elec- 
tronic (jj and electron-phonon (||| interactions. Due 
to the mainly local nature of these interactions the Dy- 
namical Mean-Field (DMF) theory seems to be an 
appropriate tool for investigations of these systems. In 
the following we will extend previous treatments |^,^|J7| 
by including both electron-phonon and electron-electron 
interactions and also by allowing for phases with long- 
range order. We will confine ourselves to the undoped 
case x = and study, in particular, the question whether 
the observed Jahn- Teller distortion is driven more by the 
electron-phonon coupling or by electronic correlations. 

Taking the limit of a very large Hund's rule coupling 
between the core 3/2 spins and the e„ electrons the re- 
sulting simplified Hamiltonian reads |8| 

H = - E *s a '(4^v + h.c.) - 1 J2 T ^ 

ijaoc' i,7 

+ \ E^-r) 2 + "14>U + 5 E &r T *r. (!) 



with 



Ti-y — E 
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c ia( c *a) creates (destructs) an electron on the ith Mn 
site in the e g orbital a =tl, with | ]) = d x 2_ y 2, | J.) = 
d 3z 2_ r 2. Because the Hund's rule coupling constant has 
assumed to be infinitely large there are no spin indices in 
H. The hopping elements t depend in the DMF theory 
on the background of the core spins. In the case of a 
ferromagnetic background t is given by the bare hopping 
element, for a paramagnetic background t contains reduc- 
tion factors H due to the double exchange mechanism. 
In the following we take for simplicity tf" — S aa itij. 
The pseudospin operator T is defined in Eq. (2) where 
the index 7 always runs only over 7 = x, z where x and z 
denote the corresponding Pauli matrices in Eq. (2). The 
effective Hubbard constant U is equal to U' — J, where 
U' is the inter-orbital Coulomb and J the exchange con- 



Jahn- Teller modes of the six oxygens surrounding the ith 
Mn site and having the local e g symmetry, dt is the time 
derivative and loq and g are the frequency and the cou- 
pling constant of the phonons. 

The Hamiltonian Eq. (1) together with the condition 
tf" = S aa 'tij, is not only invariant under the transfor- 
mations of the octahedral point group but also under 
infinitesimal rotations in orbital space. This larger sym- 
metry is due to the fact that H contains only invariants 
which are bilinear in the fields. Using the rotational sym- 
metry we can thus omit components with 7 = x in Eq. 
(1) within a good approximation, yielding Q 

H = - E ( c m c J« + H - c -) + U E n it n U 
ija. i 
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where 4n denotes the field <f>i Z . 

We want to investigate also the case where the orbitals 
order in a long-range anti-ferro orbital pattern. We there- 
fore allow for broken symmetry and divide the lattice 
into two inequivalent sublattices A and B. Within the 
DMF theory the action associated with the Hamiltonian 
(1) reduces to two single impurity actions S e ff t i where 
i = A, B and A and B stand for one site on the sublat- 
tices A and B, respectively. The two impurity actions 
read 

Seff,i= f dr f dr' Y, cUr)G^ a {r-r')c la {r') 
Jo Jo a=n 
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stant of the two e„ orbitals. 
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denote the two normal 



The two local Green's functions are related by 

GAa = Gb&, (5) 
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where a is the orbital different from a. The relations be- 
tween the Weiss field Goi a and the local Green's function 
become in the case of the Bethe lattice 



G 0Aa - lUJ v 
G 0Ba = lUJ " 



(1 — t 2 GBa 
fl — t 2 GAa- 
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Decoupling the last term in Eq. (4) by means of 
a Hubbard-Stratonovich transformation and integrating 
out the electronic degrees of freedom we obtain the new 
impurity actions 



S, 



eff,i 
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dT<p t {T)y - -j-2+ Uo)<l>i( T ) 
~ £ dr y 2 (r) + Tr ln[G$ - (gfa + VUy,)} 
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+ Tr \n[G 0l ) + (gfa + VUy,)} 
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yi{r) denotes the Hubbard-Stratonovich (HS) field. Us- 
ing Eq. (7) the local Green's functions have the exact 
representations 

G la (r -t') = J DMr)DVi(T) [G^in - r 2 ) - 

a a z a {9Un) + VUyiin^in - T 2 )]^ T ,e s °f>* /Z h (8) 

for i — A, B. Zi is the partition function. Eqs. (5)-(8) 
represent a closed system of equations for the Green's 
functions of the two coupled sublattices. 

From the equation of motion for the field 0, one ob- 
tains 

~2 (("it) _ 

(<Pb). (9) 
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Here (ni a ) is the average number of electrons in the 
orbital a on the sublattices i = A, B. The quantitiy 
(nj-f) — (riii) can ^ e considered as a staggered order pa- 
rameter. (<j)j) is the expectation value of the phonon 
field <pi. Moreover, the average number of electrons 
per manganese site is one in the undoped case, so that 
rii = (n lT +n il ) = 1. 

Assuming that the two fields <j) and y vary only slowly 
in time S e fti can be calculated by means of a gradient 
expansion |9j. The leading contribution to S e ff,i is local 
in time and is given by the following effective potential 

V,(0(r),»(r)) = ^(r) + i tf ?(r) 

- Ml " (9Mr) + VUy % {r))G 0l] {i V )] 

- e w0+ Ml + (9Ur) + Vu yi (r))G 0il (iu)]. (10) 



The lowest-order non-local contributions to S e ff t i are 
functionals of fa and yi which include also at least one 
time derivative dipi/dr or dyi/dr. 



For the evaluation of G in Eq. (8) we distinguish be- 
tween several cases. First we consider the case U = 
where the model reduces to a pure Jahn- Teller system 
and the collective coordinate y drops out. Starting from 
an educated guess for Go% a we calculate first the effective 
potentials Vi using Eq. (10) and transforming the sums 
over Matsubara frequencies into integrals along the real 
axis. Gi a , as given by Eq. (8), is calculated by using 
expansions around the extremal paths of S e ff,i- In gen- 
eral, there are two time-independent paths 4> ri f = 1,2, 
corresponding to the two minima of Vi, and one time- 
dependent instanton path. The contribution from the 
latter is very small for our considered temperatures and 
atomic masses || so we will neglect it together with the 
non-local part of S e ff t i. Around each minimum ip r we 
generate a Migdal-Eliashberg-like expansion. It is gov- 
erned by the Migdal parameter ujq jt which is rather small 
in the manganites. Detailed calculations show || that 
for all finite coupling strengths it is sufficient to include 
only the Hartree and the Fock terms. First-order vertex 
corrections to these terms turn out to be smaller by a 
factor AojoA, A being the electron-phonon coupling con- 
stant defined by A = g 2 /u> 2 t. In calculating the Fock 
term we use a phonon propagator which has been con- 
structed from the exact eigenvalues and eigenfunctions of 
the Schrodinger equation describing an atom in the effec- 
tive potential Vi . In this way anharmonic effects are fully 
taken into account. From Eq. (6) we obtain new Weiss 
fields Goi a as input for another cycle. The procedure is 
repeated until self-consistency is reached. From now on 
all energies are expressed in units of the hopping matrix 
element t, unless differently specified. 
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FIG. 1. T-A phase diagram for electrons coupled only 
to Jahn- Teller phonons. OO: ordered orbital phase; ODM 
and ODI: metallic and insulating disordered orbital phases, 
respectively. 
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Numerical results are shown in Fig. [j] where the criti- 
cal temperatures are plotted versus the electron-phonon 
coupling A. The phonon frequency has been fixed to 
luq = 0.04. Three phases can be identified: 1) an insu- 
lating phase in which the orbitals are ordered (00), 2) a 
phase in which the orbitals are disordered and the system 
is metallic (ODM), 3) a phase in which the orbitals are 
disordered and the system is insulating (OD1). The thick 
dashed line shows the metal-insulator transition temper- 
ature, the squares indicate the calculated points where 
the orbital order-disorder transition takes place. This 
is a second-order phase transition where the order pa- 
rameters {<pi) go continuously to zero as the temperature 
approaches T c from below. 

In the orbital disordered phase the two e g orbitals 
are degenerate and occupied with the same probability. 
The effective potentials Vi are symmetric with respect 
to (pi — and two regimes can be distinguished. In the 
strong coupling regime, A > A C (T), the potential shows 
two deep and well separated minima leading to a gap 
in the density of electronic states. In this case the sys- 
tem is in the insulating phase ODI. In the weak coupling 
regime, A < A C (T), the potential has only one minimum 
and the system is in the metallic phase ODM. In the or- 
dered phase the effective potential of each sublattice is 
asymmetric with respect to <pi = 0. In the weak cou- 
pling regime there is only one minimum and it is shifted 
with respect to the origin by an amount A(pi ~ {(pi)- In 
contrast to that in the strong coupling case such a shift 
is also accompanied by a lowering of one minimum with 
respect to the other minimum. The interaction with the 
phonons lifts the orbital degeneracy in an opposite way 
on the two sublattices. This creates a disproportion in 
the occupation probability of the two orbitals leading to 
the formation of a staggered pattern of occupied orbitals 
and consequently to the experimentally observed coop- 
erative distortion of the oxygen bonds. A qualitatively 
similar phase diagram has recently been given in Rcf. |1C[ | 
for the Holstein model with infinitely heavy atoms. 

Experimentally the static Jahn- Teller distortion in 
PrMn0 3 occurs around T c ~ 850-K" @ and in LaMn0 3 
around T c ~ 7hQK. Above T c these systems are insula- 
tors. If we compare these data with our phase diagram 
in Fig. [t]we obtain, taking t = 0.7eV |t3|], a value for A 
of about 5 in the insulating region. This is an extremely 
large value. In the following we show that a more realis- 
tic picture is obtained if one also includes the Coulomb 
repulsion between the e g electrons. 

Seffj depends on the two fields (pi and yt if also the 
electron-electron interaction is included. A reduction to 
one field is in this case possible if one treats the phonon 
field in the static approximation. The physical basis for 
this is the inequality loq -C t, which says that the ef- 
fective phonon mass is much larger than the electronic 
one. Performing then the following change of variables 
in Matsubara space 



y l {iuj n )VU + g<pi{0) 5 nfi = Xi{iu n )VU + A, (11) 

we can integrate out the <pi (0) coordinate and obtain the 
effective action 



\ ^2 ( 1 + JJ ~ Jf Sn '°l ^{iuj^Xii-iuJn) 



+ Tr HG-i - VTTTXxi] 
+ Tr ln[G M \ + VTT+Axj]. 



(12) 



The limit A — > correctly reproduces the action for the 
Hubbard model and the limit U —> the action for the 
static Jahn- Teller system. In the following we will con- 
centrate on the case U 3> A. Using the local approxima- 
tion for the Tr In term the effective potential reads 
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Vi( Xi (T)) = -x^rf 



T£Y" 0+ ln[l - \JTPv\xi (r)Goi-[ (iv)] 

V 

T J2e iv0+ ln[l + y/U + Xxi{T)Goii(w)]. (13) 



This means that there is no kinetic term in Eq. (12). 
Gi a can again be calculated by using expansions around 
the potential minima <p r . The propagator for the field 
x becomes instantaneous and thus can be put to zero. 
As a result only the Hartree term in each expansion 
around <p r contributes. The lowest-order corrections to 
the field propagator would come from the non-local part 
in S e tf t i and lead to a relaxational dynamics for the field 
x. Since for all couplings the damping is substantially 
smaller than lu these corrections are, however, small ||. 



0.8 



0.6 - / 
/ 



ODM 



I— 0.4 



0.2 - 



ODI 



."V 



oo 



0.0 1 — 



0.0 1.0 2.0 



3.0 4.0 

x+u 



5.0 6.0 



FIG. 2. Phase diagram showing the critical temperature 
versus X + U for a static phonon and dynamic HS field assum- 
ing U > A. 
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We have solved Eqs.(ll)-(13) in a self-consistent way. 
Figure || shows the calculated phase diagram for the case 
where the phonon, but not the Hubbard-Stratonovich 
field, is treated in the static approximation. We also as- 
sumed U A which allows to lump the two coupling con- 
stants together into an effective coupling constant A + U. 
The black dots are the calculated critical temperatures 
as a function of A + U. Our previous discussion of the 
three phases 00, ODM, and ODI in Fig. [I] also applies 
more or less to Fig. [| Remarkable is the close similar- 
ity of the phase boundaries in Figs. [I] and |^, even on a 
quantitative level. Let us estimate the Jahn- Teller cou- 
pling in the present case. Using t — 0.7eV, T c ~ O.li 
and the fact that LaMn03 is insulating above T c we ob- 
tain from Fig. || an effective coupling A + U ~ 5.3 and 
An = ua\ — hai ~ 1. The energy associated with the 
Jahn- Teller distortion can be estimated from the experi- 
mental data and is about E ~ Q.lAeV . From that, mak- 
ing use of Eq. (9) with An ~ 1, we obtain A ~ 0.3. This 
means that the experimental data can easily be repro- 
duced with U ~ i.beV and A ~ 0.3. An independent de- 
termination of U yielded 3.8eF p4| and is thus in good 
agreement with our value. The value A = 0.3 associ- 
ated with one branch of phonons is on the one hand side 
not small indicating the importance of phonons. On the 
other hand, the critical coupling strength for a metal- 
insulator transition is according to Fig. [I] about 1.5 
and thus five times larger than the obtained value. This 
means that the insulating phase above T c as well as the 
Jahn- Teller transition itself in PrMnOa and LaMnOs are 
mainly caused by electronic correlations and not by the 
electron-phonon coupling. 
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FIG. 3. Temperature dependence of the electronic den- 
sity of states on the sublattice A for the orbital a =f versus 
frequency. The left and right panels correspond to the weak- 
and strong-coupling cases, respectively. 



Fig. |3| shows the evolution of the electronic density of 
states (DOS) of one e g orbital on the sublattice A with 
temperature for two different values of the effective cou- 
pling constant A + U. The solid lines cover the tempera- 
ture range from far below the orbital ordering tempera- 
ture T c to just below T c . The dashed line shows the DOS 
just above T c . The left panel presents data in the weak- 
coupling regime A + U = 0.75. Below the critical temper- 
ature the DOS shows an asymmetry with respect to the 
Fermi level u> — 0, signalizing a disproportion of the aver- 
age electronic occupation of each e g orbital. As a result, 
ua] =/= tia[ and, at the same time, ua\ — n B [ holds. The 
difference between uai and uai increases with decreasing 
temperature. At low enough temperature a gap opens in 
the DOS at w = ReG Ba {u) +VX+U{xa)- ReG Ba {u) 
denotes the real part of (?B a (u) and (xa) the expecta- 
tion value of the field x for the action Eq. (12) on the 
sublattice A. Above T c the system is metallic. 
The right panel in Fig. || shows data for A + U = 1.8. 
Above T c the system is in the strong-coupling regime 
where the effective potential has two well separated min- 
ima causing two symmetric subbands in the DOS. With 
decreasing temperature spectral weight is transferred 
from the lower to the upper band leading to an asymme- 
try in the spectra. The gap stays practically constant un- 
til, at low enough temperatures, the orbital ordering gap 
becomes larger than the charge gap and almost all spec- 
tral weight is transferred to the right band. At T=0.04t 
we have An = ua\ — tlai ~ 0.86. 

As mentioned previously our Hamiltonian Eq. (1) is 
for the case of isotropic hopping matrix elements invari- 
ant under continuous rotations in orbital space which 
can be characterized by an ang le 6 g. As a result 9 
is not determined by Eq. (1) but by residual interactions 
such as anisotropic hoppings, higher-order invariants in 
the couplings due, for instance, to a quadratic electron- 
phonon coupling. The real situation in the manganitcs 
is even more complicated because of the observed tilt of 
the oxygen octahedra and the important role played by 
the oxygen degrees of freedom. However, we can dis- 
cuss with our Hamiltonian the temperature dependence 
of the electron density on the Mn sites if we choose, for 
instance, 9 — 7r/4 which is not far away from the experi- 
mental situation in PrMn03 and LaMnOs. Fig. ^ shows 
for A + U = 4.0 the angular dependence of the electronic 
density on the two sublattices i=A,B, defined by 

Pi (9,v) =n ir \{U0,v)\ 2 +n i l\{l\6, ( p)\ 2 , (14) 



(t \0, f) and (| \8, <p) denote the angular components 
of I T) = {d x 2_ y 2 + d 32 2_, r 2)/V2 and | |) = (d x 2_ y 2 - 
d 3z 2_ r 2)/\i / 2, respectively. 
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(15) 



(16) 



and performing the Gaussian integral over the phonon 
field we obtain 

Seff,i= f dr f dr' Y, cUr)G^ a {r-r')c la {r') 
Jo Jo a=n 



(17) 



dr VC/x i (r)[ni T (r) -n^r)] 
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FIG. 4. Angular dependence of the electronic density on 
two adjacent Mn sites for three different temperatures and 
A + U = 4.0. The vertical lobe points into the c-direction, the 
horizontal ones alternatively into the a- and fo-directions. 

The first pair of densities on the left in Fig. U cor- 
responds to the temperature T = 0.04i < T c and the 
site occupancies ua] — 0.98 and nsf = 1 — ua\ = 0.02. 
The electronic charge distribution is strongly anisotropic 
leading to alternating short and long Mn-0 bonds along 
the a- and b- directions on the two sublattices. The bond 
in the c-direction remains constant. The pair of densities 
in the middle of Fig. | corresponds to T = 0.116* < T c 
and the site occupancy n^j = 0.67. One recognizes that 
the increasing temperature decreases the anisotropy of 
the electronic distribution and thus also the distortion of 
the octahedra. Finally, the pair of densities at the right 
depicts the electron distribution on the two sublattices 
for T = 0.16t > T c and ua\ = ns] = 0.5. The electrons 
are now uniformally distributed along the three axes and 
the oxygens occupy the equilibrium positions of the cubic 
perovskite lattice. 

The above calculations for the case involving both 
the electron-electron and the electron-phonon interaction 
were based on classical, infinitely heavy atoms and the 
condition U S> A. Here we want to point out that there 
is another tractable case, that is when the fluctuations 
for the combined system are treated in a Gaussian ap- 
proximation. Performing a Hubbard-Stratonovich trans- 
formation in Eq. (4) we obtain 

Seff,i= j dr [ dr' Y < a {r)G^ a {r-T')c la {r') 

JO JO a— T J. 
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dr [V Uyi (r) - gfc (r)] [n lT (r) - mi (r)] 



where Qo = ujq y/l + X/U is the renormalized phonon fre- 
quency. We now integrate out the electronic degrees of 
freedom. The effective action depends only on x% and it 
reads 



2 ^ w 



Tr log[G^ T - VUxi] + Trlog[Go s \ + VTlx,}. (18) 



Expanding up to the first order in uj 2 we obtain 

o x ,i = --J dTXi(T) - — 

+ Trlog[G 0s 1 T -y/u + Xxi] 
+ Trlog[G 0s 1 j +VU + Xxi]. 



wg(l + U/X) dr 2 



(19) 



The action Eq. (19) is no longer a function of only the 
sum U + A but depends on both variables U and A in an 
independent way. 

The approximate action Eq. (19) reproduces cor- 
rectly the cases without phonons (A = 0) and without 
electronic correlations (U — 0). Moreover, for T ^S> 
w o\/l + U/X/2ir only the zero- frequency variable Xi(0) 
contributes to the action corresponding to the case where 
Xi can be treated as a classical and time-independent 
variable. Keeping all terms in Eq. (19) is equivalent 
to an Gaussian approximation for temporal fluctuations. 
Eq. (19) therefore represents an interpolation between 
the two exactly treated cases A = and U = which 
is valid if the temperature is not much smaller than 
the renormalized frequency u>q a/I + U/X/2n. The ac- 
tion given by Eq. (19) can be treated with the same 
methods as described above for the phonon-only case. 
Figure shows the resulting phase diagram for a fixed 
A = Ao = 0.3 and u>o = 0.04. We have chosen Ao = 0.3 as 
a representative value for the electron-phonon coupling, 
according to what we have estimated above. Compari- 
son of Figs. ^ and || indicates that treating the atoms 
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quantum-mechanically and not classically with infinite 
heavy masses yields only small changes justifying the 
above simpler treatment based on Eqs. (11)-(13). 
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Jahn- Teller transition in manganites could be made more 
realistic by including additional interactions. Anisotropic 
hopping terms, an electron-phonon interaction of second 
order, and the explicit treatment of oxygen degrees of 
freedom would remove the invariance of our Hamiltonian 
with respect to infinitesimal rotations in orbital space 
and thus determine the electronic density uniquely. Fur- 
thermore, instead of one effective electron-electron inter- 
action constant, distinct inter- and intra- orbital inter- 
action terms should be considered, which would allow to 
include also low-spin states. 
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FIG. 5. Phase diagram showing the critical temperature 
versus U + Xq, with Ao = 0.3, for dynamic phonons and HS 
fields in a Gaussian approximation. 

In conclusion, we have extended previous treatments of 
manganites within the DMF theory including also phases 
with long-range order and, in addition to the electron- 
phonon coupling, the electron-electron interaction. Us- 
ing expansions around stationary paths of the action 
which are controlled by the Migdal parameter we were 
able to treat the weak- and the strong-coupling cases on 
an equal footing. However, since the temporal fluctua- 
tions of the phonon and the Hubbard-Stratonovich (HS) 
fields are different we considered three special cases: a) 
the phonon-only case; b) the case where the phonons are 
treated statically, the HS field dynamically; and c) treat- 
ing both fields dynamically in the Gaussian approxima- 
tion. For each case the phase diagram was shown. Apply- 
ing the results to undoped manganites it was concluded 
that the observed Jahn- Teller transition is mainly driven 
by the electron-electron and only to a lesser extent by the 
electron-phonon interaction. Our simple model for the 
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